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We refute an often invoked theorem which claims that a periodic orbit with an odd number of 
real Floquet multipliers greater than unity can never be stabilized by time-delayed feedback control 
in the form proposed by Pyragas. Using a generic normal form, we demonstrate that the unstable 
periodic orbit generated by a subcritical Hopf bifurcation, which has a single real unstable Floquet 
multiplier, can in fact be stabilized. We derive explicit analytical conditions for the control matrix 
in terms of the amplitude and the phase of the feedback control gain, and present a numerical 
example. Our results are of relevance for a wide range of systems in physics, chemistry, technology, 
and life sciences, where subcritical Hopf bifurcations occur. 
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The stabilization of unstable and chaotic systems is a 
central issue in applied nonlinear science 0j S ■ Start- 
ing with the work of Ott, Grebogi and Yorke j^, a vari- 
ety of methods have been developed in order to stabilize 
unstable periodic orbits (UPOs) embedded in a chaotic 
attractor by employing tiny control forces. A particularly 
simple and efficient scheme is time-delayed feedback as 
suggested by Pyragas [1] . It is an attempt to stabilize pe- 
riodic orbits of minimal period T by a feedback control 
which involves a time delay r = nT , for suitable positive 
integer n. A linear feedback example is 



z{t) = .f{\-4t)) + B[z{t-T)-z{t)] 



(1) 



where z{t) = /(A, z{t)) describes a d-dimensional nonlin- 
ear dynamical system with bifurcation parameter A and 
an unstable orbit of period T . B is & suitably chosen 
constant feedback control matrix. Typical choices arc 
multiples of the identity or of rotations, or matrices of 
low rank. More general nonlinear feedbacks are conceiv- 
able, of course. The main point, however, is that the 
Pyragas choice rp = nT of the delay time eliminates 
the feedback term in case of successful stabilization and 
thus recovers the original T-periodic solution z{i). In this 
sense the method is noninvasive. Although time delayed 
feedback control has been widely used with great success 
in real world problems in physics, chemistry, biology, and 
medicine, e.g. d, 0, i, i, M llH, E, E, Q III, El, O [ll , 
severe limitations are imposed by the common belief that 
certain orbits cannot be stabilized for any strength of the 
control force. In fact, it has been contended that periodic 
orbits with an odd number of real Floquet multipliers 
greater than unity cannot be stabilized by the Pyragas 
method [11, [H, |2l|, [H, [H, Q , even if the simple scheme 
([T]) is extended by multiple delays in form of an infinite 
series (25l |. To circumvent this restriction other, more 
complicated, control schemes, like an oscillating feedback 
[2^ . or the introduction of an additional, unstable degree 
of freedom [24,[27|, have been proposed. In this letter, we 
claim, and show by example, that the general limitation 



for orbits with an odd number of real unstable Floquet 
multipliers greater than unity does not hold, but that 
stabilization may be possible for suitable choices of B. 
We illustrate this with an example which consists of an 
unstable periodic orbit generated by a subcritical Hopf 
bifurcation, refuting the theorem in [20| . 

Consider the normal form of a subcritical Hopf bifur- 
cation, extended by a time delayed feedback term 

z{t) = [A + j + (1 + Z7)k(t)n z{t) + b[z{t - r) - z{t)]{2) 

with z G C and real parameters A and 7. Here the 
Hopf frequency is normalized to unity. The feedback 
matrix B is represented by multiplication with a com- 
plex number b = ba + ibj = boe^^ with real bii,bi,(3, 
and positive bo. Note that the nonlinearity f{X,z{t)) = 
[X + i + (1 + ij)\z{t)\'^'\ z{t) commutes with complex ro- 
tations. Hence the Hopf bifurcations from the trivial so- 
lution z = at simple imaginary eigenvalue rj = iuj ^ 
produce rotating wave solutions z{t) = z(0) exp ii'^t) 
with period T even in the nonlinear case and with delay 
terms. This follows from uniqueness of the emanating 
Hopf branches. 

Transforming Eq. ^ to amplitude and phase variables 
r, 9 using z{t) = r(t)e*^'*-', wc obtain at 6 = 



r{t) = (A + r^)r 
e{t) = l + 7r2. 



(3) 
(4) 



An unstable periodic orbit (UPO) with = —A and 
period T = 27r/(l - 7A) exists for A < 0. At A = a 
subcritical Hopf bifurcation occurs. The Pyragas control 
method chooses delays as rp = nT. This defines the local 
Pyragas curve in the (A, T)-plane for any n 6 N 



TP (A) 



27rn 
1 - 7A 



27rn(l + 7A + . . . ) 



(5) 



which emanates from the Hopf bifurcation point A = 0. 
Under further nondegeneracy conditions, the Hopf point 



2 



A = 0; r = nT (n G No) continues to a Hopf bifurcation 
curve T/f (A) for A < 0. We determine this HopJ curve 
next. It is characterized by purely imaginary eigenvalues 
77 = iuj of the transcendental characteristic equation 



r] = X + i + b (e-''^ - 1) 



(6) 



which results from the linearization at the steady state 
z = of the delayed system 

Separating Eq. ^ into real and imaginary parts 

= A + 6o[cos(/?-cjt) -cos/?] (7) 
w-l = 6o[sin(/3 - t^r) - sin/3] (8) 

and using trigonometric identities to eliminate 
yields an explicit expression for the multivalued Hopf 
curve th{\) for given control amplitude 60 and phase 
/3: 



± arccos 



(9) 



1-60 sin (3 ± A / A(26o cos /3 - A) + 6^ sin^ [3 



Note that th is not defined in case of /? = and A < 0. 
Thus complex 6 is a necessary condition for the existence 
of the Hopf curve in the subcritical regime A < 0. Fig. [1] 
displays the family of Hopf curves, n G No, Eq. ([9]), and 
the Pyragas curve n = 1, Eq. ([5|), in the (A, r) plane. In 
Fig. [ijb) the domains of instability of the trivial steady 
state 2 = 0, bounded by the Hopf curves, are marked 
by light grey shading (yellow online). The dimensions 
of the unstable manifold of z = are given in parenthe- 
ses along the r-axis in Fig. [ijb). By construction, the 
period of the bifurcating periodic orbits becomes equal 
to Tp ~ nT along the Pyragas curve, since the time- 
delayed feedback term vanishes. Standard exchange of 
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FIG. 1: (Color online) Pyragas (red dashed) and Hopf (black 
solid) curves in the (A, r)-plane: (a) Hopf bifurcation curves 
n = Q, 10, (b) Hopf bifurcation curves 71 = 0, 1 in an en- 
larged scale. Yellow shading marks the domains of unstable 
z = and numbers in parentheses denote the dimension of the 
unstable manifold of z = (7 = —10, 60 — 0.3 and /3 — 7r/4). 

stability results [2^ , which hold verbatim for delay equa- 
tions, then assert that the bifurcating branch of periodic 
solutions locally inherits linear asymptotic (in)stability 
from the trivial steady state, i.e., it consists of stable 
periodic orbits on the Pyragas curve Tp(A) inside the 



yellow domains for small |A|. Note that an unstable triv- 
ial steady state is not a sufficient condition for stabi- 
lization of the subcritical orbit, but other (e.g., global) 
bifurcations at A < must be considered as well. More 
precisely, for small |A| the unstable periodic orbits pos- 
sess a single Floquet multiplier /i = exp(AT) G (l,oo), 
near unity, which is simple. All other nontrivial Floquet 
multipliers lie strictly inside the complex unit circle. In 
particular, the (strong) unstable dimension of these peri- 
odic orbits is odd, here 1, and their unstable manifold is 
two-dimensional. This is shown in Fig. [21 which depicts 
solutions A of the characteristic equation of the periodic 
solution on the Pyragas curve. Panel (a) (top) shows the 
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FIG. 2: (Color online) (a) top: Real part of Floquet exponents 
A of the periodic orbit vs. feedback amplitude 60. bottom: 
Real part of eigenvalue r; of steady state vs. feedback am- 
plitude bo. (b): Floquet multipliers /i — exp(Ar) (red) in 
the complex plane with the feedback amplitude bo G [0, 0.3] 
as a parameter, (c): radii of periodic orbits. Solid (dashed) 
lines correspond to stable (unstable) orbits. (A — —0.005, 



7 = -10, T^jt 



13 = Tv/4). 



dependence of the real part of the critical Floquet expo- 
nent A on the amplitude of the feedback gain bo. The 
largest real part is positive for bo = 0. Thus the periodic 
orbit is unstable. As the amplitude of the feedback gain 
increases, the largest real part of the eigenvalue becomes 
smaller and eventually changes sign. Hence the periodic 
orbit is stabilized. Note that an infinite number of Flo- 
quet exponents are created by the control scheme; their 
real parts tend to —00 in the limit 60 0, and some 
of them may cross over to positive real parts for larger 
60 (blue curve), terminating the stability of the periodic 
orbit. Panel (b) of Fig. [2] shows the behavior of the Flo- 
quet multipliers fi = exp(AT) in the complex plane with 
the increasing amplitude of the feedback gain bo as a pa- 
rameter (marked by arrows). There is an isolated real 
multiplier crossing the unit circle at /i = 1, in contrast to 
the result stated in [2^ . This is caused by a transcritical 
bifurcation (TC) in which the subcritical Pyragas orbit 
(whose radius is given by r = (— A)^/^ independently of 
the control amplitude bo) collides with a delay-induced 
periodic orbit, as shown in Fig. [2l[c). This delay-induced 
orbit is generated at a finite value of the control am- 
plitude bo (SN) by a saddle-node bifurcation (collision 
with another unstable delay-induced periodic orbit). At 
TC, the subcritical orbit and the delay-induced orbit ex- 
change stability. The latter vanishes at a subcritical Hopf 
(subH) bifurcation at which the trivial steady state be- 
comes unstable. Except at TC, the delay-induced orbit 



3 



has a period T ^ t. Note that for small 69 the subcritical 
orbit is unstable, while z = is stable, but the respective 
exchanges of stability occur at slightly different values of 
60, corresponding to TC and subH. This is also corrob- 
orated by Fig. [IJa) (bottom), which displays the largest 
real part of the eigenvalues 77 of the steady state z = 0. 
The possible existence of such delay-induced periodic or- 
bits with T ^ which results in a Floquet multiplier 
/i = 1 of multiplicity two at TC, was overlooked in [20|. 

Next we analyse the conditions under which stabiliza- 
tion of the subcritical periodic orbit is possible. From 
Fig. [Hb) it is evident that the Pyragas curve must lie in- 
side the yellow region, i.e., the Pyragas and Hopf curves 
emanating from the point (A, r) = (0, 27r) must locally 
satisfy the inequahty th{X) < Tp(A) for A < 0. More 
generally, let us investigate the eigenvalue crossings of 
the Hopf eigenvalues rj = ioj along the r-axis of Fig. [1] 
In particular we derive conditions for the unstable dimen- 
sions of the trivial steady state near the Hopf bifurcation 
point A = in our model equation ([2|). On the r-axis 
(A = 0), the characteristic equation ([6]) for 77 = iuj is 
reduced to 



rj = i + b 



-1 



and we obtain two series of Hopf points given by 



0<t' 



< T, 
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27rn 

2/3 + 27rn 
1 - 2&0 sin/3 



(n = 0,l,2,...). 



(10) 

(11) 
(12) 



The corresponding Hopf frequencies are u ~ 1 and 
uj^ = 1 — 26q sin (3, respectively. Note that series A con- 
sists of all Pyragas points, since = nT = In the 
series B the integers n have to be chosen such that the 
delay r^f > 0. The case bo sin/3 = 1/2, only, corresponds 
to lli^ — and does not occur for finite delays t. 

We evaluate the crossing directions of the critical Hopf 
eigenvalues next, along the positive r-axis and for both 
series. Abbreviating -^rj by r]r the crossing direction is 
given by sign(Re i]r)- Implicit differentiation of (fTO| with 
respect to r at 77 = implies 

sign(Re rjr) ~ — sign(a;) sign(sin(wr — /3)). (13) 

We are interested specifically in the Pyragas-Hopf points 
of series A (marked by red dots in Fig.l) where r = 
= 2-Kn and lo ~ uo^ = 1. Indeed sign(Re 77,-) = 
sign(sin/3) > holds, provided we assume < /3 < tt, 
i.e., bi > for the feedback gain. This condition alone, 
however, is not sufficient to guarantee stability of the 
steady state for r < 2777r. We also have to consider the 
crossing direction sign(Re 77r) along series B, uj^ = 1 — 
2&osin/3, Lo^T^ = 2P + 2Tm, for < /3 < tt. Eq. ^ now 
implies sign(Re rjr) = sign((2&o sin/3 — l)sin/3). 

To compensate for the destabilization of z = upon 
each crossing of any point r^ = 2^71, wc must re- 
quire stabilization (sign(Re 1]^-) < 0) at each point 



of series B. This requires < /3 < arcsin(l/(26o)) or 
TT — arcsin (l/(2&o)) < /3 < tt. The distance between two 
successive points r^ and r,f_^i is 2ti /lo^ > 2n. Therefore, 
there is at most one r^ between any two successive Hopf 
points of series A. Stabilization requires exactly one such 
r^, specifically: r^_-^ < t^_i < for all fc 1, 2, . . . ,n. 
This condition is satisfied if, and only if, 



< /3 < /3: 



(14) 



where < /3* < tt is the unique solution of the transcen- 
dental equation 



TT 



2nbo sin/3* ~ 1. 



(15) 



B 

fc-1 



This holds because the condition t^_i < r^ 
fails when t^_i = t^. Eg. pi]) represents a necessary but 
not sufficient condition that the Pyragas choice rp = nT 
for the delay time will stabilize the periodic orbit. 

To evaluate the second condition, th < Tp near 
(A, r) = (0, 27r), we expand the exponential in the charac- 
teristic eq. ^ for lot w 2Tm, and obtain the approximate 
Hopf curve for small |A|: 

th{X) ^2nn- ^{2TrnbR + l)X. (16) 



Recalling the Pyragas stabilization condition 

th(A) < rp(A) is therefore satisfied for A < if, and 
only if. 



2^) < 



(17) 



Ea. p?|) defines a domain in the plane of the complex 
feedback gain b = bp + ibj = b^e^^ bounded from below 
(for 7 < < &/) by the straight line 



^ —7 27r77,) 
Eq. P3|) represents a curve 60 (/3), i.e.. 



6n- 



277 sin /3 \ TT 



(18) 



(19) 



which forms the upper boundary of a domain given by the 
inequality (Hi]). Thus (HH) and (dH) describe the bound- 
aries of the domain of control in the complex plane of the 
feedback gain b in the limit of small A. Figl3] depicts this 
domain of control for 77 = 1, i.e., a time-delay r = j-z^- 



The lower and upper solid curves correspond to Eq. (|18[) 
and Eq. respectively. The color code displays the 

numerical result of the largest real part, wherever < 0, 
of the Floquet exponent, calculated from linearization of 
the amplitude and phase equations around the periodic 
orbit. Outside the color shaded areas the periodic or- 
bit is not stabilized. With increasing |A| the domain of 
stabilization shrinks, as the deviations from the linear 
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approximation (jl6p become larger. For sufficiently large 
A| stabilization is no longer possible, in agreement with 
FigHfb). Note that for real values of 6, i.e., /3 = 0, no 
stabilization occurs at all. Hence, stabilization fails if the 
feedback matrix i? is a multiple of the identity matrix. 



mcinschaft in the framework of Sfb 555. 




Re(b) 



FIG. 3: (Color online) Domain of control in the plane of the 
complex feedback gain b = boe''^ for three different values of 
the bifurcation parameter A. The black solid curves indicate 
the boundary of stability in the limit A / 0, see (|18fl . (|19fl . 
The color-shading shows the magnitude of the largest (nega- 
tive) real part of the Floquet exponents of the periodic orbit 



(7 = -10, 
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In conclusion, we have refuted a theorem which claims 
that a periodic orbit with an odd number of real Flo- 
quet multipliers greater than unity can never be stabi- 
lized by time-delayed feedback control. For this purpose 
we have analysed the generic example of the normal form 
of a subcritical Hopf bifurcation, which is paradigmatic 
for a large class of nonlinear systems. We have worked 
out explicit analytical conditions for stabilization of the 
periodic orbit generated by a subcritical Hopf bifurca- 
tion in terms of the amplitude and the phase of the feed- 
back control gain [2^ . Our results underline the crucial 
role of a non-vanishing phase of the control signal for 
stabilization of periodic orbits violating the odd number 
limitation. The feedback phase is readily accessible and 
can be adjusted, for instance, in laser systems, where 
subcritical Hopf bifurcation scenarios are abundant and 
Pyragas control can be realized via coupling to an exter- 
nal Fabry-Perot resonator The importance of the 
feedback phase for the stabilization of steady states in 
lasers [l^l and neural systems js^l, as well as for stabi- 
lization of periodic orbits by a time-delayed feedback con- 
trol scheme using spatio-temporal filtering [sij , has been 
noted recently. Here, we have shown that the odd num- 
ber limitation does not hold in general, which opens up 
fundamental questions as well as a wide range of applica- 
tions. The result will not only be important for practical 
applications in physical sciences, technology, and life sci- 
ences, where one might often desire to stabilize periodic 
orbits with an odd number of positive Floquet exponents, 
but also for tracking of unstable orbits and bifurcation 
analysis using time-delayed feedback control ^3^ . 
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